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ABSTRACT 

This  technical  memorandum  describes  a  mathematical  model  for 
determining  the  expected  number  of  men  suppressed  by  a  rapid  fire 
machine  gun  mounted  on  a  helicopter.  Also  included  is  a  description 
of  the  computer  program  of  the  model  for  the  IBM  7030  (STRETCH). 
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FOREWORD 

The  mathematical  model  and  computer  program  described  in  this 
technical  memorandum  were  developed  as  part  of  WEPTASK  RM3773-109/ 
210-1/F008-99-08. 

The  programming  was  done  by  Miss  J.  A.  Knight  under  the  super¬ 
vision  of  Mr.  R.  0.  Brancolini. 
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INTRODUCTION 


The  model  described  in  this  report  determines  within  a  specified 
target  area  the  expected  number  of  men  suppressed  by  a  rapid  fire 
machine  gun  mounted  on  a  helicopter.  This  model  is  intended  to  evaluate 
the  effectiveness  of  various  armament  systems  in  suppressing  a  number  of 
personnel  targets. 

For  the  purpose  of  this  model,  suppression  is  defined  as  the  expect¬ 
ation  of  impacting  a  predetermined  number  of  projectiles  (or  more)  within 
a  predetermined  rectangular  or  square  area  around  a  man.  Both  the  pre¬ 
determined  number  of  bullets  for  suppression  and  the  size  of  the 
rectangle  or  square  can  be  varied  in  the  model  presented  below,  i.e., 
the  definition  of  suppression  can  be  varied. 

This  model  has  been  programmed  for  the  IBM  7030  (STRETCH)  and  a 
listing  of  the  computer  program  is  presented  in  Appendix  E.  Besides 
the  program  listing,  the  input  format  and  sample  of  the  output  are 
presented  in  Appendix  B  and  Appendix  C,  respectively. 

The  problem  for  which  the  model  was  developed  and  the  assumptions 
upon  which  the  model  is  based  are  presented  in  the  following  sections. 

PROBLEM  STATEMENT 


The  problem  for  which  the  model  was  developed  can  be  stated  as 
follows:  Given  the  number  of  rounds  per  burst,  the  dispersion  (down 
range  and  cross  range)  of  the  mean  point  of  the  impacts  of  the  machine 
gun  bursts,  the  dispersion  (down  range  and  cross  range)  of  the  bullets 
in  a  burst,  and  the  number  of  men  in  the  target  area,  determine  the 
expected  number  of  men  suppressed  in  a  flat  target  area  of  specified 
size. 

ASSUMPTIONS 


The  following  assumptions  are  used  in  the  model: 

1.  The  target  area  is  flat  and  the  terrain  cover  is  of  such  a 
nature  so  as  not  to  deflect  the  bullets. 

i 

2.  The  men  are  randomly  dispersed  throughout  the  target  area. 

3.  The  areas  of  suppression  of  the  men  do  not  intersect. 

4.  The  target  area  is  rectangular  or  square. 

5.  The  line  of  fire  of  the  machine  gun  is  always  parallel  to  one 
of  the  sides  of  the  target  area. 
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6.  The  mean  points  of  the  impacts  (MPIs)  of  the  bursts  have  an 
independent  bivariate  normal  distribution. 

7.  The  bullets  within  a  burst  form  an  independent  bivariate 
normal  distribution. 

8.  The  down  range  and  cross  range  dispersions  of  both  the 
bullets  within  a  burst  and  the  MPIs  of  the  bursts  are  known. 

9.  The  number  of  bullets  per  burst  is  constant  and  known. 

MODEL  DESCRIPTION 


Although  the  following  description  of  the  model  is  presented  in 
only  enough  detail  to  make  it  understandable,  an  effort  has  been  made 
to  keep  the  mathematical  symbols  of  this  presentation  consistent  with 
the  symbols  used  in  the  computer  program  coding. 

First  the  target  area  is  divided  into  rectangular  or  square  blocks 
whose  dimensions  are  those  of  the  area  which  defines  suppression.  (See 
hatched  area  of  Figure  1.)  These  dimensions  are  Ax  in  the  direction 
perpendicular  to  the  line  of  fire  and  Ay  in  the  direction  parallel  to 
the  line  of  fire. 

Since  it  is  possible  for  the  mean  point  of  impact  (MPI)  of  a  burst 
to  lie  outside  the  target  area  and  still  have  bullets  from  that  burst 
strike  within  the  target  area,  it  may  be  necessary  to  add  extra  blocks 
of  the  same  size  as  the  target  blocks  around  the  target  area.  (See 
Figure  1  unlined  blocks.)  These  additional  impact  blocks  are  added  until 
the  probability  of  having  an  MPI  in  one  is  small  (on  the  order  of  .0005 
or  less).  The  number  of  blocks  to  be  added  (if  any)  is  dependent  upon 
the  dispersion  of  both  the  bursts  (MPIs)  and  the  bullets. 

The  Impact  area  consists  of  these  additional  blocks  plus  the  target 
blocks.  If  the  dispersions  of  the  MPIs  are  such  that  the  probability 
of  getting  an  MPI  outside  the  target  area  is  less  than  .0005,  the  target 
area  and  the  Impact  area  are  identical.  The  impact  area  is  always 
square  or  rectangular,  and  the  number  of  additional  impact  blocks  (if 
any)  in  the  x-  and  y-directions  are  predetermined  and  are  inputs  to  the 
computer  program. 

Next,  labeling  of  the  blocks  begins  with  the  lower  lefthand  corner. 
This  first  block  is  labeled  (1,1).  The  first  coordinate  indicates  the 
ith 

row,  and  the  second  coordinate  indicates  the  column  of  the  impact 
area.  All  impact  blocks  are  labeled  in  this  way.  The  target  blocks 
are  labeled  by  the  coordinate  pair  of  the  impact  block  they  represent 
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*Target  Area  and  Impact  Area 


*  The  hatched  area  is  the  target  area,  and  the  blocks  within  it  are 
target  blocks. 

The  total  area  --  lined  and  clear  --  is  the  impact  area,  and  the 
blocks  within  it  (including  the  lined  blocks)  are  impact  blocks. 


Figure  1. 
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(the  Impact  area  ia  always  at  least  as  large  as  the  target  area).  To 
avoid  confusion  when  referring  to  target  blocks,  k  will  designate  the 
column  of  the  target  block,  and  1  will  designate  the  row. 


Let 

■  standard  deviation  of  the  MPI  of  the  weapon  burst  cross 
range 

‘’Vw  ■  standard  deviation  of  the  MPI  of  the  weapon  burst  down 
range 

It  is  assumed  that  a  cartesian  coordinate  system  with  its  origin  at 
the  centroid  of  the  impact  area  is  overlayed  on  the  impact  area.  The 
x-direction  is  cross  range;  the  y-direction  is  down  range.  Using  this 
coordinate  system,  let 


nix  -  mean  of  the  distribution  of  MPIs  cross  range 

my  -  mean  of  the  distribution  of  MPIs  down  range 

P(i,J)  ■  probability  that  the  MPI  of  a  given  burst  is  in  Impact 
block  (i,j). 

Then: 


P(i,j) 


1 

*yv 


-1/2 

e 


dxdy 


since  x  and  y  are  assumed  independent. 


|u-v|  -Ax 

| t-s j  -Ay,  i.e.,  u  and  v  are  the  cross  range  and  s  and  t  are  the 
down  range  distances  from  the  mean  to  the  side  of  impact 
block  (i, j) . 

The  method  of  determining  the  values  of  u,  v,  s,  and  t  adds  nothing 
to  the  presentation  of  the  model,  and  it  is  not  included.  However,  the 
method  can  be  found  in  the  program  listing  of  Appendix  E. 
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Let 


■  cross  range  standard  deviation  of  the  bullets  for  any 
5  HP  I 

-  down  range  standard  deviation  of  the  bullets  for  any 
^  MPI 

The  mean  of  the  distribution  of  bullets  is  assumed  to  be  at  the  centroid 
of  the  impact  block  for  each  impact  block.  That  is,  the  mean  of  the 
bullets  is  assumed  to  be  (0,0)  in  a  cartesian  coordinate  system  whose 
origin  is  the  center  of  the  impact  block. 

Let 


Pik,l)/(i,j)  -  probability  that  a  bullet  hits  a  target  block 
L  (k,l)  given  an  impact  in  impact  block  (i,j) 


Then 


1 

2lf0ry%  *yb 


''d  rg 

Jc  Jf 


2 


"I/2  <#- 

xb 


dxdy 


M  -Ax 

td- c|  ■  Ay  >  l.e.,  g  and  £  are  the  cross  range  and  c  and  d  are 
the  down  range  distances  from  the  center  of  impact 
block  (i,j)  to  the  sides  of  target  block  (k,l). 

Again  the  method  of  determining  g,  f,  and  d,c  adds  nothing  to  the 
presentation  and  is  not  given.  It  too  can  be  determined  from  the  program 
listing  of  Appendix  E. 

Let 

N  =  the  number  of  bullets  per  weapon  burst 

B  l(k, 1)  / (i , j)  |  =  expected  number  of  bullets  in  target  block 
J  (k,l)  given  the  MPI  in  impact  block  (i,j) 

B  [(k,l)/(i,j)J  =  N-P  [(k,l)/(i,j)j 
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Since  suppression  can  occur  only  in  those  blocks  receiving  an 
expected  number  of  bullets  equal  to  or  greater  than  some  tolerance  q, 
it  is  necessary  to  count  and  record  for  each  impact  block  (i,j)  the 
number  of  target  blocks  (k,l)  for  which  B  f(k,  1)  /(i, j)l  2?  q.  For  each 
impact  block  (i,j)  call  this  quantity  T(i,j). 

T(i,j)  •  number  of  blocks  in  which  suppression  can  occur  given  an 
impact  in  impact  block  (i,j) 

K  -  number  of  targets  in  the  target  area 

B  *  number  of  target  blocks  into  which  the  target  area  is  divided 

£ 

”  =  probability  that  a  target  is  in  any  given  target  block 

B 

E(i,j)  “  expected  number  of  targets  suppressed  by  an  MPI  in 
impact  block  (i,j) 

E(i,j)  =  S  •  T(i,j) 

D 

E  =  expected  number  of  targets  suppressed  by  a  single  machine  gun 


E  =  E  E  /e  ( i ,  j )  •  P(i,j>) 

all  j  all  i  V  J 

=  E  E  (f  '  T(i,j)  *  P(i,jA 

all  j  all  i  ) 


=  |  •  E  E  T(i, j)  •  P(i, 


all  j  all  i 


j) 


In  the  case  where  the  tolerance  q=l  the  following  option  can  be 

used : 


Determine  the  expected  number  of  bullets  striking  the  target  area 
but  not  accounted  for  in  the  target  blocks  which  form  T(i,j)  for  each 
impact  block  (i,j).  Call  this  expected  number  of  bullets  w'(i,j). 

w'  (i,  j)  =  Eb  ^(k,  1)  /i ,  j  )  J  ,  for  all  B  £(k,l)/i,  j)  jcq 

Let 


w(i,j)  =  integer  portion  of  w'(i,j) 
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Let 


-  T(i,j)  +  w(i,j) 


E*(i,j)  ■  expected  number  of  targets  suppressed  by  an  MPI  in 
impact  block  (i,j)  using  the  option 

E'd.j)  -  |  •  T’(i,j) 

C*  ■  expected  number  of  targets  suppressed  using  the 

option  by  a  single  machine  gun  burst 

E'  -  |  •  £  £  T'd.j)  -  P(i,j) 

all  j  all  i 

This  option  is  reasonably  accurate  only  when  q  ■  1.  For  other  values 
of  q  it  is  not  usable. 

COMPUTER  PROGRAM 

The  mathematical  model  described  in  the  previous  section  has  been 
programmed  for  the  IBM  7030  (STRETCH).  The  programming  language  is 
FORTRAN  IV. 

The  symbols  used  in  the  computer  program  are  defined  in  Appendix  A, 
and  the  inputs  to  the  program  are  indicated.  Appendix  B  presents  a 
sample  input  sheet;  Appendix  C  gives  an  example  of  the  output  from  the 
program.  Appendix  D  is  a  flow  chart  of  the  computer  program,  and 
Appendix  E  is  a  listing  of  the  program. 

Some  of  the  input  data  such  as  angle  (angle  of  descent),  velocity, 
altitude,  and  range  from  target  are  not  actually  used  in  the  computer 
program,  but  they  are  included  in  the  output  for  identification  in  those 
cases  where  these  data  are  used  to  compute  the  various  dispersions.  If 
for  a  particular  theoretical  study  these  data  are  not  applicable,  zeroes 
can  be  used  in  the  input. 

RESULTS  AVAILABLE  FROM  THE  MODEL 

This  section  describes  how  the  model  can  be  used. 

The  parameters  that  can  vary  are: 

(1)  ^xw  ^yw  =  cross  range  and  down  range  dispersions, 

respectively,  of  each  weapon  burst 


11 


(2)  d'xb  an(*  =  cross  range  and  down  range  dispersions, 

respectively,  of  each  weapon  burst 

(3)  Ax  and  Ay  ■  dimensions  of  target  and  impact  blocks  and 

also  area  of  suppressions 

(4)  q  ■  expected  number  of  bullets  needed  in  target  block  to 

suppress  a  man 

(5)  N  *  number  of  bullets  per  weapon  burst 

A  change  in  any  of  these  parameters  will  produce  a  change  in  E,  the 
expected  number  of  men  suppressed. 

A  study  can  be  performed  on  the  effect  of  changing  any  two  parameters 
by  determining  E  for  each  value  of  the  changing  parameters  and  holding  all 
other  parameters  constant.  Figure  2  shows  a  possible  family  of  curves 
for  changes  in  and  0^  .  (These  curves  are  presented  only  for  the 

purpose  of  illustration,  ana  their  form  has  no  relation  to  those  an  actual 
study  would  produce.) 
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Illustration  of  Parametric  Study 


(°xb  and  <7^  vary) 


Ax 


Ay 


REMAIN 
I  CONSTANT 


VARIOUS 

'  b 


It  is  possible  to  get  a  family  of  such  families  by  giving  those 
parameters  that  remain  constant  new  values. 


Figure  2. 


13 


DEFINITION  OF  SYMBOLS  USED  IN  COMPUTER  PROGRAM 


*  I  RUN 

*  NHAA 

*  NHVK 

*  NHAF 

*  NHRF 

*  NTX 

*  NTY 

*  NAX 

*  NAY 

*  XDEL 

*  YDEL 

*  T0L 

*  N 

*  K 

*  XMI 

*  YMI 

*  SIGIX 

*  SIGIY 

*  SIGBX 

*  SIGBY 

*  I0P 


=  Run  Indicator  J  =  pother  run  to  follow 

1  =  Last  run 

=  Helicopter's  Angle  of  Descent 
=  Helicopter's  Velocity 
=  Helicopter's  Altitude 
=  Helicopter's  Range  from  Target 

*  Number  of  Target  Cells  in  x-direction 
=  Number  of  Target  Cells  in  y-direction 

=  Number  of  additional  cells  in  x-direction 
=  Number  of  additional  cells  in  y-direction 
=  Cell  dimension  in  x-direction 

*  Cell  dimension  in  y-dit"ci:ion 

*  Expected  number  of  bullets  needed  in  each  cell  for 
suppression 

=  Number  of  bullets  in  each  machine  gun  burst 
=  Number  of  targets 

=  Mean  of  MPIs  of  distribution  in  x-direction 
=  Mean  of  MPIs  of  distribution  in  y-direction 
=  Standard  deviation  of  MPIs  of  distribution  in  x-direction 

=  Standard  deviation  of  MBI^  of  distribution : in  y-idirection 

=  Standard  deviation  of  bullets  of  distribution  in  x-direction 

*  Standard  deviation  of  bullets  of  distribution  in  y-direction 

=  Option  ?  =  option 

1  =  Include  option 


*  These  quantities  are  program  inputs. 


A-l 


NRUN 
NIX 
NIY 
LMAX 
KMAX 
XTAR 
YTAR 
XIMP 
YIMP 
CONST  1 
2 

3 

4 

**  5 
NAXP1 

NAYP1 

TOLDN 

HX 

IIY 

*ERF(T) 


Run  number 

Number  of  Impact  cells  In  x-direction 
Number  of  Impact  cells  In  y-direction 
Index  number  of  last  target  block  in  y-direction 
Index  number  of  last  target  block  in  x-direction 
x-dimension  of  target  area 
y- dimension  of  target  area 
x-dimension  of  impact  area 
y-dlmension  of  impact  area 
\ 


V^Wx 

v*s 


Constants  needed  in  transforming  the 
normal  distribution  to  ERF* 


-  S 

B 

■  Index  number  in  x-direction  of  target  block  in  lower  left- 
hand  corner 

*  Index  number  in  y-direction  of  target  block  in  lower  left- 
hand  corner 

*  Cutoff  point  in  terms  of  probabilities 

*  Distance  in  x-direction  from  (0,0)  of  coordinate  system  at 
center  of  impact  area  to  the  edge  of  the  impact  area 

=  Distance  in  y-direction  from  (0,0)  of  coordinate  system  at 
center  of  impact  area  to  the  edge  of  the  impact  area 


=  JL 


2 

~z  dz 


**See  text,  Model  Description,  for  definition  of  symbols. 
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BX  ■  Number  of  blocks  in  x-dlrection  excluding  the  center  block 

for  the  bullets  in  a  burst.  If  3X  is  not  an  integer,  then 
IX  makes  it  an  integer. 

BY  ■  Number  of  blocks  in  y-direction  excluding  the  center  block 

for  the  bullets  in  a  burst.  If  BY  is  not  an  integer,  then 
IY  makes  it  an  integer. 

IXP1  *»  Number  of  blocks  in  x-direction  including  the  center  block 
for  the  bullets  in  a  burst. 

IYP1  ■  Number  of  blocks  in  y-direction  including  the  center  block 
for  the  bullets  in  a  burst. 

PX  =  Probability  that  an  MPI  occurs  along  a  given  line  segment 

in  the  x-direction. 

PY  ■  Probability  that  an  MPI  occurs  along  a  given  line  segment 

in  the  y-direction. 

PXP  *  Probability  that  a  bullet  strikes  along  a  given  line  segment 

in  the  x-direction. 

PYP  ■  Probability  that  a  bullet  strikes  along  a  given  line  segment 

in  the  y-direction. 

T(I,J)  *  Number  of  target  blocks  which  meet  suppression  requirement 
given  an  impact  in  impact  block  (I,J) 

TP(I,J)  =  Number  of  target  blocks  for  the  option  which  meets  the 

suppression  requirement  given  an  impact  in  impact  block  (I,J) 

IW  =*  Integer  portion  of  those  bullets  expected  to  strike  the  target 

area  in  target  blocks  which  do  uot  meet  the  suppression  require 
ment. 

E  *  Expected  number  of  targets  suppressed 

EP  -  For  option  expected  number  of  targets  suppressed 
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L'n 


DATA  CARD  LAYOUT 


HE1IC0PTER  ARMAMP'T  STUDY 


RUN  NUMflFR  1 

ANGLE  -  45  CEGREES  ALTITUDE  =  90GC  FELT 

VELOCITY  »  120  KNOTS  RANGE  *  10CC  TE6T 

TARGET  AREA  IS  40. CC  FEET  ev  42. CO  Ft F T 
NUMBER  CF  TARGETS  IS  14 

TARGET  DIMENSIONS  ARE  10. CC  FEET  BY  t . 2C  FEFT 


RFAN  OF  MEAN  POINTS  OF  IMPACT  IN  X-0  H-  CCT  TUN  IS 
MEAN  OF  MEAN  POINTS  OF  IMPACT  IN  V-GTRfcTION  IS 


-P.  FtET 
11.  FFFT 


STANDARD  DEVIATION  OF  MEAN  POINTS  GF 
STANDARC  DEVIATION  CF  MEAN  POINTS  CF 

STANDARD  DEVIATION  CF  BULLETS  IN  EACr 
RTANOAKD  DEVIATION  OF  PULLFTS  IN  EACH 


IMPACT 

IN 

X— C! RECT ION 

IS 

10.  -u  0 

EFF  1 

IMPACT 

IN 

V-7IRECTICN 

IS 

«  •  JOG 

FFFT 

-<1>R  ST 

IN 

X-RIRFCTICN 

IS 

e.cce 

FEET 

RIJR  ST 

IN 

V— DIRECTION 

IS 

b.'JCC 

FEPT 

NURriro  CF  PULLETS  PER  BURST  IS  4b 
rxPCTED  NUMHFR  OF  SUPPRESSIONS  IS  ■  3.7 

CPTll.N 

L*P  (.Trn  N JPR I R  OF  SUPPRESSIONS  IS  *  4 .► 


Uf.  .  I  f.r,  i  l  in  SECONDS  * 


.757 


I 


,  — —  .  ■ - 

Compute  the  probability  that 
a  bullet  hits  a  target  block 
for  all  target  blocks  given 
an  impact  in  the  first,  second, 
etc . ,  target  block 


For  each  target^\^ 
block  of  the  given  impact 
block  is  the  probability  : 
\±he  required  tolerancg^ 


Compute  the  expected  number  of 
bullets  striking  the  target 
area  but  not  accounted  for  in  the 
target  blocks  >.  the  tolerance 


Sum  of  total  number  of 
target  blocks  £  the 
required  tolerance 


For  the  given  impact  celL 
compute  the  expected  number 
of  targets  suppressed 


Have  all 

impact  cells  been 
^v^cons  i  de  red^^-^ 
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TYPE « COMP 1LGO* FORTRAN, PM 


SU0TYPE»FOPTRAN,LMAP*PUNCH 
DIMENSION  T(  2*500)*  TP(2500) 

COMMON  I  PUN *  NHA A  * NHVK , NHAF  , NHRF , NTX *  NT y  «  NAX  » NAY * XDEL « YDEL  «  TOL  « 

1  N , K  «  XM l  «  YM I  ,SIGIX,SIGI Y , S I GPX « S I GBY »  IOP«NRUN 

NPiJNsO 

1  READ  2,  I  RUN*  NW AA  , NHVK  «  NHflc , NHRF ,  NTY  * NT Y  «  NAX  *  NA  Y  ♦ XDEL  *YDEL»TOL« 

1  N*K«XMI  ♦ YMI  , E I G I  X. S IG I Y«  S IGBX«S IGBYt  IOP 

nrijn=nrun+ 1 

3  FORMAT  (  I 2«  3Xt I 3«  3X«  I *«  ?X •  I  *  «  2X. 14 . 3X« ! 3  «2X« I 3«  2X«  I  3«2X« I  3  «  2X« 

1  F5,2*2X,F5.2,3y*F6.3*3X, 13/ 

2  I 3,2X,F»0.0*3X,F1 0 . 0  * 2X , F 7 . 3 , 2X . F7 . 3 . 2X « F7 . 3 , 2X , F7 . 3 , 

3  3X«  I  1  »3X*F4.2 ) 

CALL  SET  IT 

NIX*  NTX+?*NAX 

N  I  V  =  NTY+2*NAY 

I- ( (NI X*NI Y ) .GT.250C >  STOP 

CALL  HAE(T.TP,NIX.NIY) 

I  r  (  I  RIJN.LT  •  1  )  CO  TO  1 

RE  T>jRN 

END 

SUBTYPE* FOPTRAN.LMAP.^UNCH 
SUEROUT  INF  (T  ,  T»  ,MX  ,NI  v  ) 

r I MFNSI ON  T (NI X  *NI Y )  , TP ( M  y ,NI V  > ,PX (5FD )  ,PY (50  0  )  «PXP (500 ) ♦ PYP ( 50  0  ) 
COMMON  I  PUN, NHA A , NHVK , NHA-  , NHRF ,Ntv,NtY.NAX,NAY«XDEL*YDEL»TOL* 

1  N.  <  «  XM  I  ,  YM  I  ,  F  I  G  I  X  ,  S  I  G  I  Y  ,  S  1  Oi'X  «  S  J  0 n Y  *  I  OP  «  NRUN 


LMAX*NTY  +  NAY 


<MAX*NTX  +  NAX 

XT AP  *  XDEL*FLOAT<NTX) 

yjAR  =  YDEL*FLOAT(NTY) 

XIMP  =  XTAR+2 • *FLOAT (NAX )*XDEL 
YIMP  =  YTAR+2.*FL0AT (NAY )*YDEL 

fn=n 

C0N5T1 =1 .4 1 421 35623731 *SIG IX 
C0NST2=1 .4142135623731 4SIGI Y 
C0NRT3=1  .41421  3 e 62373  1  *S  I  GBX 
C0N3T4=1  .4  14213?  S2373 1*51  GBY 
rON5T^= (FLOAT (K  5  )/ (FLOAT (NTX*NTY )  ) 

NAXP1 =NAX+1 
tolon=tol/float (N ) 

HX=X I MP/2. 

HY=v I MP/2. 

RX= ( 3.5*S IGRX ) /X0EL-.5 
RY= ( 3 • 5*S I GBY ) /YDEL- • 5 
=  NJ  A  V  4-  1 

°P  I  NT  2  1  «  NRUN.NHAA  ,  NH  A  F  «  NHV<  ,  NHR~  ,  x'TAP»yTAP#K«  XDEL  «  YDEL  «  XM  I  iYM| 
21  fodvat  (  i  H2 , 26X,  25HHFL  I  COPTFR  ARMAMENT  r  T'JIY / / 

1  1 OHOPuN  NtjMOFQ  ,11/ 

3  IHO.fiy.BHANGLF  =  P.PH  PFGREFS.  1  4X«  1  1  HAL.T  I  TUDE  =  I  4  «  5H  FEET/ 

3  *>Y.  1  1  Hvf-'LOr  I  TY  r  (4,  AH  KNOT r- «  1  ?X  « RHPANOF  =  I  4  «  fh  FEr  T  // 

4  1H0«6X* 1 FHTAPGFT  APEA  IC  FR,2,QH  F r E T  EY  F8.P.5H  FEET/ 

j  7X.21HNUMRFP  OF  T  APGPT E  I c  13/ 
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7X«22HTARGET  DIMENSIONS  ARE  F5.2*9H  FEET  BY  F5.2*5H  FEET// 


7  1  HO »  6X « 48HMEAN  OF  MEAN  POINTS  OF  IMPACT  IN  X-DIRECTION  IS  FlO*0* 

8  5H  FEET/7X  «44HMEAN  OF  MEAN  POINTS  OF  IMPACT  IN  Y-D I RECT I  ON « 

9  AH  IS  F  1  0  •  0  «  5H  FEET//) 

PRINT  22 »  S I G I  X i S I G I Y  »  S I GBX *  S  I GB Y  «  N 
22  FORMAT ( 1  HO «  6X «  4 6H STANDARD  DEVIATION  OF  MEAN  POINTS  OF  IMPACT  IN, 

1  16H  X-DIRECTION  IS  F7.3*5H  FEET/ 

2  7X « 46HST  ANOARD  DEVIATION  OF  MEAN  POINTS  OF  IMPACT  IN* 

3  16H  Y-D I RECT I  ON  IS  F7.3«5H  FEET// 

4  1  HO , 6X , 46H STANDARD  DEVIATION  OF  BULLFTS  IN  EACH  BURST  IN* 

5  16H  X-DIRECTION  IS  F7.3,FH  FEET/7X* 13HST ANDARD  DEVIATION* 

6  43H  OF  BULLETS  IN  EACH  BURST  IN  Y-D 1 RECT I  ON  IS  FB.3.5H  FEET// 

7  1  HO  «  6X  «  3 1 HNUM9ER  OF  BULLETS  PER  BURST  jc  13) 

X=— HX-XMI 

PI =.5*EPF (ABS (X/CONS1  1  )  ) 

DO  F  I = 1  *  N I  X 
V=  Y-t-XHFL 

P?r  ,«*Fpe-  (  flq<;  (x/CON^T  1  )) 

!cOO  4,4,2 
7  ir(v-vDFL)  3,3*4 

3  PX (  I  )  =P 1 +P2 
PI  =n? 

GC  TO  F 

4  PX  (  I  )  =  ARS  (  P  1  -P”1  ) 

PI  =P2 

F  CONTINUF 
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Ys-HY-YMI 


PI*  «5#ERF ( ABS ( V /CONST? ) ) 
D010  I*1«NIV 

y*y+ydel 

P 2=  .5*ERF(ABS(Y/C0N?T2 )) 
IF(Y)  9  «  9  « 7 

7  IF(Y-YDEL)  8«8«9 

8  PY ( I ) =P 1 +P2 
PI  =P£ 

GO  TO  10 

9  PY (  I  )  =ABS (PI -P2 ) 

PI  =P2 

10  CONTINUE 

IX  =RX+,9Q9 
IY  =  B  Y+  • 999 
Y  I  =  I  Y 

VI  =  IX 

nx  =- ( X  I  +  #5 ) *XDEL 
nv  ~  -  (y  i  +  »c'  j*ydel 

X  =  DX 

L  =  ?* IX  +  2 

PI = .F*ERF ( ABP ( X/CCNPTT  )  ) 
IXPl  = I X+ 1 

no i r  i * i « ixpi 

X  =  x-*-y  PPL 

PP-.^^FPFlAP^fX/CONPT1)  ) 
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PXP< l )  *  ABS<P1*P2) 


PI  *P2 
L«L-1 

PXP (L )*PXP( 1 ) 

15  continue 

PXP< IXP1 1 *2##P2 
V  »dy 

L*2*lY+2 

PI «.5*ERF(ABS  ( Y/CONST 4 )  ) 

I YP 1 « I Y+ 1 
D020  I  =  1  «  I YP 1 
Y* Y+YOEL 

P2  « . 5*ERF ( ABS ( Y/C0NST4 ) ) 

P YP ( I  )=ABS<P1-P2> 

PI  =  P2 
L  =  L-1 

PYP(L)=PYP( I ) 

20  CONTINUE 

pYP ( 1 YP1  ) *2  «  *P2 
0060  J=1 «NI  Y 
D060  1  =  1  *  N I  X 
T ( I , J)=  0 
=  0 

TP  (  I  «  J  )  =0 

DO  50  L=NAYP1 «LM AX 

IF  (  (  IA9S(L-J)-I YPl  )  .GE.O  )  GO  TO  25 
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LI  *  IYP1  +  IABS< J-L  ) 

P2  »PYP(L1  ) 

GO  TO  30 
25  P2  *  0 

30  DO  50  K=NAXP1 *KMAX 

IF ( ( IABS<K-I)-IXP1 ) .GE.O )  GO  TO  35 
LI  *  I XP1  +  I ABS (K-I  ) 

PI «PXP (LI  ) 

GO  TO  40 
35  P 1 =0 
40  n  =Pl*P2 

IF (P.LT.TOLDN)  GO  TO  45 
T ( I « J)  =T( I « J)+l • 

GO  TO  50 
45  5UM  =  SUM+P 
"  CONTINUE" 

\'H-  FUM*FN 

=  I  W 

TP (I « J) *T ( I ♦ J )+W 

A '  FONT  I NUc 
r  -  o 

DO  6*  J= 1  «N I  V 
P?=°Y ( J ) 
no  55  I=1 ,NIX 
° 1 =PX (  I  ) 

C  =  r  +  I  ,  J)*P1  -4  =  ? 
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65  CONTINUE 

F«  E*C0NST5 
PRINT  66 «E 

I F ( IOP.EO.0  )  GO  TO  75 

66  FORMAT  (37H0EXPE  *TED  NUMBER  OF  SUPPRESSIONS  IS  *  F6*l> 

EP  =  0 

no  70  j=i ,niv 

P2  =  PY ( J  ) 

DO  70  1*1  » N  I  X 
PI =PX ( I  ) 

frp  =  FP+TP(  I  «  J)#P1#P2 
70  CONTINUE 

EP=EP*C0NST5 
PRINT  71 ,EP 

7J  FORMAT  (7H0OPTIGN/  07H0FXPECTED  NUMRFR  OF  SUPPRESSIONS  IS 
75  CALL  INTVL(TFFC) 

PRINT  BO.TSEC 

FORMAT  (P7HORUNNINC-  TIME  IN  SECONDS  =  »F7.3) 

PFTl  ION 
r'ND 

T  SyRTYPE  *  DATA 
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